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METHOD OF SHIM OPTIMIZATION 

BACKGROUND OF THE INVENTION 
Field of Invention 

5 This invention relates to a method of optimizing the effects of a shim on magnetic field 

in an imaging apparatus or a multi-sample spectrometer; and more particularly, to such a method 
for determining values of currents applied to the coils of the shim so as to enable homogeneous 
magnetic field to be applied to a plurality of cells of a probe used therein. 
Discussion of the Related Art 

10 In the art, there is provided an imaging system, such as those utilizing nuclear magnetic 

resonance (NMR), wherein a radio frequency (RF) probe having one or more cells therein is 
used to obtain data from one or more samples. Each cell has a space for holding a sample, such 
as a mouse. The plurality of cells are disposed within a space of a magnet, which could be a 
superconducting magnet, an electro magnet or a permanent magnet, which provides a static 

15 magnetic field within which the cells are disposed. The RF signal is provided, for example, by a 
common RF power source to each cell, through a power divider, which transmits an RF signal to 
each cell. The resulting imaging signal is received by each cell and is coupled through a 
transmit/receive (T/R) circuit to a separate receiver, detector and data storage memory, all of 
which are known in the art. A plurality of probes, each having one or more cells, can also be 

20 used. 

The static magnetic field, contained within the sample space, usually requires adjustment 
for desired properties, such as magnetic homogeneity. This can be done by a shim system of 
one or more coils which are located at different positions with respect to the magnetic field to 
produce magnetic field gradients used to compensate for spatial variations of the magnetic field. 

25 The shim coils need to be adjusted to produce uniform magnetic field value Bo to all of 

the probe cells in order to produce optimal results. 

One method for shim adjustment is to combine the RF output signals for each cell and 
couple the combined RF signals to a common receiver. With use of a simulated sample, such as 
a simulated mouse, in each cell, the plurality of shim coils can be adjusted to obtain overall 

30 signal height or narrowest overall line width, as desired. Another method, and similar to the 

foregoing method, is to combine the signals after detection and use the resulting signal height or 
line width to determine the optimal shim setting. In order to avoid problems with signal 
phasing, the absolute signal amplitude can be used instead of a phase adjusted absorption mode 
response. 
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In liquid spectroscopy, wherein multiple probe cells are used with multiple shim coils for 
adjustment of common signals, similar problems exist. 

The foregoing shim adjustment methods have all been found to have deficiencies and 
leave much to be desired. In these cases, adjustment requires repeated optimization steps. Also, 
5 adjustment of one shim coil changes the optimization required of other shim coils so that it is 
difficult to obtain the desired degree of optimization for all of the shim coils without requiring 
an undue number of readjustment steps to achieve the desired optimization. Also, there is 
continued uncertainty of whether the best fit or adjustment is ever obtained. Thus, there is, in 
the art, uncertainty when an adjustment procedure should be terminated. 

10 

SUMMARY OF THE INVENTION 

Accordingly, an object of the invention is to overcome the aforementioned and other 
deficiencies and disadvantages of the prior art. 

Another object is to optimize the effects of the shim system on magnetic fields to provide 
15 uniformity of magnetic field values at different probe cells disposed therein. 

A further object is to obtain the same value of resonant frequency of simulated samples 
at different probe cells disposed in a static magnetic field. 

A further object is to provide optimal shim settings in a liquid spectroscopy system 
utilizing a probe having a plurality of cells. 
20 Another object is to optimize shim settings so that a minimum amount of errors occurs in 

a magnetic field at different probe cells. 

The foregoing and other objects are attained in the invention method which is based on 
the mathematical relation of the frequency in a cell Fj depends on the linear relationship with 
currents l u h, In of the shim coils 1, 2, ... N, with a constant term Gj, the frequency being with 
25 no currents in the shim coils. This linear relationship is expressed by the following: 

N 

F j= G j + Z a jkV 

k=l 

One object of the invention is to find the shim current values that will result in the same 
resonant frequency at each cell. 

This is done by first measuring the resonant frequency of the cells when no current is 
30 applied to the shim coils giving the values of Gj, and then, for each cell j measuring the change 
of the resonant frequency (Fj - Gj) when a current is applied to one shim coil k without any 
current being applied to the other shim coils. The change in frequency (Fj - Gj) divided by the 
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current Ik that produced that change gives a sensitivity factor ajk of the one cell of the selected 
shim coil k. For the selected cell j the process is carried out for the current in each shim coil k. 
The process is repeated for each cell j giving a matrix of sensitivity factors ajk. The shim 
currents values required to produce the same frequency in each cell are calculated by a matrix 
5 multiplication of the inverse of the sensitivity matrix ajk and difference of the target frequency 
and initial frequency of each cell Fj - Gj. 

For the case, when the number of RF probe cells M is greater than the number of shim 
coils N, the matrix of values is operated on to obtain the least square of errors thereby 
minimizing the difference of the magnetic field values from the desired values Fj - Gj. 

10 A feature of the invention is a method of optimizing effects of shims on adjusting a 

magnetic field applied to a plurality of cells of a probe, comprising the steps of determining a 
value of current which is applied to a plurality of shim coils so that the resonant frequency of 
each cell of the probe is substantially the same and so that the plurality of cells of the probe are 
subjected to a uniform magnetic field. 

15 Another feature is the method wherein the value of the current is determined using the 

following steps: measuring resonant frequency of each of the probe cells disposed in the 
magnetic field without any current being supplied to the plurality of coils of the shim, thereby to 
obtain Gj, G 2 , G M , wherein M is the number of the probe cell; measuring resonant frequency 
of a first cell of the probe using a current Ii applied to a first shim coil and without any current 

20 being applied to the other shim coils to obtain Fh, wherein the first subnumber is the number of 
the probe cell and the second subnumber is the number of the shim coil which is energized by 
current Ii; measuring resonant frequency of a second cell with current Ii being applied to a first 
shim coil and without any current being applied to the other shim coils, to obtain F21. This 
process is contained for the remaining cells up through Fmi- The process is then repeated 

25 measuring the frequencies of the cells for a current h in shim coil 2, with no other shim currents 
present. Continuing this process through all shim coils leads to a matrix of values Fjk. The 
sensitivity matrix ajk is then found using the following: 

■j*— — 

wherein the matrix becomes: 



3 
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A = 



a ll a l2 
a 2l a 22 

a Ml a M2 



a 2N 
a MN 



The frequencies at each cell j, Fj and current in each shim k, Ik in matrix form are: 
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5 With these definitions the linear relationship between the frequencies of the cells, F, and 

the shim currents, I, is expressed in matrix form by the following: 

F=G+AI 

10 One object of the invention is to find the shim current values, which produce the same 

frequency, or the best approximation to the same frequency, in each cell. 

A feature is the method of optimization applied when the number of independent shim 
coils N is equal to the number of probe cells M, wherein the matrix is square and has a unique 
inverse A' 1 . To find the shim current values I k , which produce the same frequency in each cell 

15 the following matrix equation is solved by setting all values of Fj equal, i.e. set Fi = F 2 = . . . = 
Fm- 

I = A-' (F-G), 

20 Wherein A" 1 is the inverse matrix of A, i.e. A" 1 A =1 . 

Another feature is the foregoing method being utilized when the number of shim coils N 
is equal to or greater than the number of probe cells M. 

A further feature is the optimization method when the number of shim coils N is less 
than the number of probe cells M, and wherein the currents in each shim are calculated 
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according to a least square mathematical operation, and the least error in the magnetic field 
values are obtained for each cell. 

Another feature is wherein the method is used to optimize shims used in adjustment of 
magnetic fields in an imaging system wherein the magnetic field is supplied by a magnet, and 
5 the plurality of probe cells are disposed in the magnetic field and the shims are disposed close to 
the cells and within the magnetic field. 

Another feature is use of the method in a liquid spectroscopy apparatus wherein the 
probe cells and the shim coils are utilized for a plurality of samples. 

10 BRIEF DESCRIPTION OF THE DRAWINGS 

For a better understanding of the present invention, reference is made to the 
accompanying drawings: 

Fig. 1 shows a schematic view of assembly comprising a magnet with shim coils 
disposed therein, shim coils comprise a plurality of RF probe cells, each having a sample 
15 therein. 

Fig. 2 shows a schematic view of assembly comprising a magnet with shim coils 
disposed therein. Shim coils have a plurality of RF probe cells, wherein each RF probe cell has 
a plurality of samples. 

20 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

This invention is applicable to an imaging system, such as a nuclear magnetic resonance 
imaging system, wherein one or more probes, each having one or more cells therein, are 
disposed in a static magnetic field, which may be supplied by a superconductive magnet, an 
electromagnet, or a permanent magnet. To adjust or account for various inhomogeneities, such 

25 as residual magnetism, imperfections in the magnetic coils, etc, of the static magnetic field, one 
or more shims, each having one or more coils, such as linear and quadratic gradient and higher 
order shim coils, are utilized within the magnetic field, as known in the art. Further, the shim 
coils are made to produce magnetic field corrections that are linearly independent, i.e., the 
corrections produced by anyone of the shim coils cannot be achieved by some combination of 

30 currents in the remaining shim coils. As above discussed, there are various methods known in 
the art, for optimizing the effects of the shim coils on the magnetic field, but, none of which are 
entirely satisfactory. 

The instant invention intends to optimize the setting of the shim currents in a reliable and 
simple manner. 
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The invention is also applicable to high resolution liquid spectroscopy systems, wherein 
a plurality of the cells are utilized concurrently to analyze a plurality of samples. 

It was discovered that for a probe having a plurality of cells (each usable for a separate 
sample, e.g. a mouse or a simulated sample, e.g. a simulated mouse) and disposed in a static 
5 magnetic field with a shim system having a plurality of shim coils to adjust the static magnetic 
field to compensate for various inhomogeneities, the best overall Bo magnetic field homogeneity 
for all the probe cells is obtained when the same value of magnetic field is obtained at each cell. 
This was obtained when the resonant frequency of each cell, e.g. with a simulated sample 
therein, had the same value. 
10 The resonant frequency Fj of a particular probe cell j is equal to the resonant frequency 

Gj of that probe cell due to the static magnetic field and without any influence due to current 
being applied to any of the shim coils, plus the sum of the contributions to the resonant 
frequency made by current applied to the other shim coils, that is: 

F j= G j + i>jkV 

15 The following description describes a method of determining shim coil settings that give 

the best overall measurement to minimize errors in the magnetic field values of the cells. 

When using a probe having a plurality of cells M, and a shim system having a plurality 
of coils N within the magnetic field, a first step in the invention method is to measure the 
resonant frequency of each cell without any current being applied to any of the shim coils. The 

20 result can be called Gi, G2, Gm, wherein M is the number of the cell being measured for 
resonant frequency. 

Next, the resonant frequency of each cell Fj is measured with a current Ii being applied to 
shim coil 1 with no current being applied to the other shim coils. The current Ij can be of any 
non-zero arbitrary value. The resonant frequency measurement is repeated for each cell j 
25 yielding Fji. The following equation is attained: 

F jl - G 3 



Then, the resonant frequency of the first cell is measured (called Fi 2 ) with a current I2 of 
non-zero arbitrary same value applied to shim coil 2, with no current being applied to the other 
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shim coils. Doing this for all cells produces a factor a j2 from the above equation (2), wherein the 
subscript 1 is replaced by the subscript 2. 

The foregoing steps are repeated for each of the probe cells and for each of the shim coils 
using currents I k so as to produce a plurality of factors ajk, for each row M (equal to number of 
5 cells) and each column N (equal to number of coils), which then produces a matrix product A, as 
below: 



10 



A = 



a ll a 12 

a 21 a 22 



a lN 
a 2N 



a Ml a M2 



.(3) 



The matrix form for frequencies of the cells Fj, the residual frequencies Gj, and currents 
in each shim coil Ik are as set forth below: 



F = 







/c > 




h 




G = 




I = 















Equation (1), written in matrix form is: 

F-G + AI ...(4) 

15 wherein F is the resonant frequency of all cells, G is the resonant frequency of all cells in the 

static magnetic field without current being applied to the coils, and A is the matrix product of 

the factors matrix (3) due to influence of current in the shim coils. 

When N = M, i.e. the number of shim coils is equal to the number of probe cells, the 

matrix A is square and has an inverse matrix A* 1 . As mentioned above the shim coils must be 

20 linearly independent, with the result that no row of matrix A is a linear combination of other 

rows. This condition is satisfied if the determinant of A does not vanish. If the determinant of 

A = 0, i.e. | A | =0, then the coefficients of ay are not linearly independent, and the matrix A 

does not have an inverse. This can happen if two of the shim coils produce identical fields at 

each cell and therefore are redundant. When the number of non-redundant shim coils are equal 

7 
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10 



15 



to the number of probe cell the matrix A has an inverse. The inverse matrix A -1 satisfies the 
condition A" 1 A =U, wherein U is the unit matrix with unity on the diagonal and O's elsewhere. 

The calculation of the inverse of a matrix is well known in the art and does not form part 
of this disclosure. Methods of analyzing matrices are found in many books, for example 
"Matrix Analysis" by Roger A. Horn and Charles R. Johnson. Various software programs such 
as Mathcad by MathSoft, Inc. Cambridge, Massachusetts are available that numerically invert a 
matrix and perform other matrix operations. 

When the number of shim coils is less than the number of probe cells there is insufficient 
number of controls to guarantee the same field can be created at each cell. In this case, the least 
square operation is used to come as close as possible. 

The problem is to find the shim currents I required to produce the same field Fo at each 
cell, i.e Fj = F 2 , = F M = F 0 . The currents I are then obtained by the following equation: 

I-A-'CF-G) ...(5) 

The current is thus a matrix of currents obtained in equation 5 by using the following 
frequency matrix: 

I = . with F = • ••• (6) 



12 



When N is greater than M, there are more shim coils than are required to produce 

identical resonant frequencies in each cell. 
20 In that case, one can select M shim coils only for use in the shimming of the apparatus. 

Usually, lower order shim coils values are used. In that case, the foregoing rules for the case of 

N being equal to M, would also be applicable. 

Other different numbers for N coils can be used, and the results suitably utilized. For 

example, data from all of the coils can be utilized to optimize some other parameters, such as 
25 total shim power, or to improve homogeneity of one of the cells. Another procedure is to place 

more than one sample in some or all of the cells so the number of samples is less than or equal to 

the number of shim coils. In this case the number of shim coils would be greater than the 

number of probe cells, but need not be greater than the number of samples. 
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10 



When N is less than M, the number of shim coils N is insufficient to guarantee that the 
frequencies of all M probe cells would be of the same value. Thus, it was found that the shim 
coil settings could be obtained by using the least square of errors of frequencies from the desired 
frequency Fo. That is, a method of least squares is used to obtain a current Ik that best matches 
the frequency of each cell Fj to desired frequency F 0 . In such a method, the squares of the errors 
at each cell are added up to generate an error function, namely: 



M / * M ( N V 



(7) 



To find the N current values that minimize E, we set the derivatives dE/dlj, dE/dl2, . dE/dlu 
equal to zero. For the case of dE/dI p , this becomes: 

...(8) 



dE/ 



The right hand equality can be expressed in the following matrix form: 

A T AI = A T (F - G) 



(9) 



15 



Here, A T is the transposition of matrix A, I is a 1 by N the column matrix, and F - G is a 
1 by M column matrix, as follows: 



F - G = 



F 0 " Gl 
F 0 - G 2 



F 0 " G M 



(10) 



The values of A and G were determined by measurements described above. Fo is the desired 

frequency of operation. Since F 0 can be selected to be of any desired value, it is assumed that 
20 current I applied to one of the shim coils actually controls the main magnetic field, i.e. , it 

applies the same magnetic field to all cells of the probe. 

To determine the value of shim currents which provide the best fit, equation (9) must be 

solved for the shim current I. This is done by matrix multiplication of equation (9) by (A T A) _i . 

This yields the following equation: 
25 I = (A T A) _I A T (F - G) ...(11 ) 
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It was assumed that N is less than M. In this case, neither A nor A is a square; however, 
A T A is a square and has an inverse (A T A) _1 . If A were square, then an inverse A" 1 would exist, 
and equation (11) would be reduced to equation (5), which is the case for when N is equal to M. 

Advantageously, with the invention a simple set of measurements provides results which 
5 are used to provide a current values which can be applied to the set of shim coils to produce the 
same resonant frequencies of all the probe cells when N = M and N is greater than M. Also, in 
another embodiment when N is less than M, these measurements are operated on using the least 
square method to give the current values that produce the least error of the frequencies for all the 
cells, thus enabling homogeneous magnetic field to be applied to all the cells. In the foregoing 
10 manner, the invention provides a simple and reliable method for providing uniform and 
homogeneous magnetic field to all cells of a probe. 

The foregoing description is illustrative of the principles of the invention. Numerous 
modifications and extensions thereof would be apparent to the worker skilled in the art. All such 
modifications are to be considered to be within the spirit and scope of the invention. 
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